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In this paper we revisit the issue of determining the oscillating primordial scalar power spectrum 
and oscillating equation of state of dark energy from the astronomical observations. By performing 
a global analysis with the Markov Chain Monte Carlo method, we find that the current observations 
from five-year WMAP and SDSS-LRG matter power spectrum, as well as the "union" supernovae 
' sample, constrain the oscillating index of primordial spectrum and oscillating equation of state of 

, dark energy with the amplitude less than |namp| < 0.116 and |Wamp| < 0.232 at 95% confidence level, 

. respectively. This result shows that the oscillatory structures on the primordial scalar spectrum and 

the equation of state of dark energy are still allowed by the current data. Furthermore, we point 
out that these kinds of modulation effects will be detectable (or gotten a stronger constraint) in the 
near future astronomical observations, such as the PLANCK satellite, LAMOST telescope and the 
currently ongoing supernovae projects SNLS. 

I. INTRODUCTION 

o, . . 

^ Recent advances in observational cosmology have revealed that our universe has experienced at least two different 
stages of accelerated expansion. One is the inflation in the very early universe when its tiny patch was superluminally 
stretched to become our observable Universe today. This can naturally explain why the universe is flat, homogeneous 
I and isotropic. Inflation is driven by a potential energy of a scalar (or multi-scalar) called inflaton and its quantum 
O fluctuations turn out to be the primordial density fluctuations which seed the observed large-scale structures (LSS) 
^ , and anisotropy of cosmic microwave background radiation (CMB). The other one is accelerating expansion driven by 
^ ' dark energy (DE) which dominates the energy density of the universe currently. Understanding the nature of dark 
energy is among the biggest problems in modern physics and has been studied actively in the literature. 

At present, the high quality observational data, CMBP, 0, [l], LSS[1] and type la supernovae (SN Ia)[| and so 
on, have provided the stringent constraints on cosmological parameters. For example, current data can constrain the 
primordial scalar power spectrum index Us and the energy density of dark energy component f2de to 1% level f6|. 
Besides the current observations, there are many ongoing projects, such as PLANCK[3], LAMOStQ and SNLS[9|. 
These projects will provide more accurate measurements on CMB temperature anisotropics and polarization, LSS 
p\| ' matter power spectrum and the luminosity distance, which will be helpful for studies of inflation and dark energy and 
' \ determinations of cosmological parameters. 
' Generally, the current data analysis bases on the simple parameterizations of the primordial power spectrum and 
the equation of state (EoS) of DE. However, we note that some inflation models can generate the power spectrum 
. with some modulated wiggles. This picture can be realized by the inflaton fleld with a step function of the potential 
^ ' [13, [HI or oscillating potential [l^]- The effects from the Trans-Planckian initial condition can also lead to oscillations 
• ^ „ which do imprint directly on the primordial scalar power spectrum[l3, [13) [Hi- The bouncing model driven by the 
Quintom[l^ matter can also give rise to some wiggles in the primordial scale-invariant power spectrum, because in this 
^ model the universe initially experiences a contracting stage, after the contracting phase it bounces to an inflationary 
■ - - phase, therefore the primordial fluctuations in sub-hubble region would deviate from that generated in Bunch-Davies 
vacuum [l7j. The Quintom bounce provide a solution to initial singularity problem, and in this scenario there's no 
Trans-Planckian problem since we can choose the initial condition via contracting phase. The featured primordial 
scalar perturbation spectrum with local bumps are studied in Ref.[18]. 

In some sense, the study for DE is similar to inflation, either in model building or data fltting. The featured EoS of 
DE, especially the oscillatory behavior EoS can provide us an unconventional evolution of universe. Observationally, 
these kinds of modulated EoS will leave clew on the bubble diagram or the matter power spectrum as well as the 
temperature power spectrum of CMB, which give us some hints to test such a scenario. The periodic oscillatory 
EoS can be realized by the two scalar field Quintom matter [Toj. And in Ref . [20| . a class of Quintom models with an 
oscillating equation of state have been studied. In such a scenario, the early inflation and the current acceleration 
of the universe can be unified, the scale factor keeps increasing from one period to another and leads naturally to a 
highly flat universe. The periodic recurs will not lead to big crunch nor big rip and the coincidence problem can be 
reconciled by the vicissitudinary repetition. Also, the studied relevant to these kind of oscillating DE can be found 
in [111, [H, [11 [11. 
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Given the above progress in the astronomical observations and physical motivations, it can potentially lead us to 
study issues related to the featured structure of primordial spectrum and the EoS of DE. In this paper, we aim to 
study the constraints on such kind of featured parametrization of the primordial spectrum as well as the EoS of DE. 
We perform the global fitting by using the Markov Chain Monte Carlo method with the current data from CMB, LSS 
and SN, and also from the simulated future projects, such as PLANCK, 5-year Supernovae Legacy survey (SNLS) 
and LAMOST. 

The paper is organized as follows: In section II we introduce the method of fitting and data we used in our analysis. 
The results and discussions are given in section III. The last section is on the summary. 



II. METHOD AND DATA 



As mentioned before, there are many theoretical models which can provide the non-trivial structure on the primordial 
power spectrum and EoS of DE. In order to study the oscillating primordial scalar power spectrum independent of 
the specific model, we parameterize the power spectrum Py as: 



lnPyik) = lnA,{k) + [nso{ko) 
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where, is the amplitude of primordial power spectrum, ko is the scale pivot and is set as 0.05 Mpc h 



characterizes the tilt of spectrum while namp and n{„ 
the traditional definition of spectral index, we get 



denote the contribution of featured oscillation. Comparing with 
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and one can find that for small rttrc, this parametrization goes back to the traditional form, i.e. the last term in the 
right hand side of equation above gives rise to the running of the scalar spectral index. 
For the parametrization of DE, we take the form given in 25,] : 



w{a) =Wo + Wamp siu [wfro ln(a)] , 



(3) 



This oscillating behavior in the EoS can lead to the modulations on the Hubble diagram or a recurrent universe which 
unifies the early inflation and the current acceleration. In Refs. (2ll . [2^ . [2^ [27j some preliminary studies have been 
presented on this kind of DE model. We noticed that there are some hints on the oscillating behavior for example 
as shown in the Fig. 10 of the SNIa paper [ll] and see also paper in Ref. [1^. Our sine function has the advantage 
of preserving the oscillating feature of the DE EoS at high redshift measured by the CMB data. For simplicity and 
focusing on the study at lower redshift, we set Wfro to be 37r/2 in order to allow the EoS to evolve more than one 
period within the redshift range of z = to z = 2, where the SNIa data are most robust. For the ACDM model. 
Wo = -1 and Wamp = 0. 

We extend the publicly available MCMC package CosmoMC^ 30]by including dark energy perturbation[3l|, and we 
assume the adiabatic initial condition in our calculation. The most general parameter space is: 



(Wh , , Bs , T, Wo , Wamp , "-sO , J^amp , ^frc , ) , 



(4) 



where w;, = Q,i,h? and ujc = Vlch? with f2f, and Q,c being the physical baryon and cold dark matter densities relative 
to the critical density, Og is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at 

A, 



^amp ; 



decoupling, t is the optical depth to re- ionization, wo and Wamp are the EoS parameters of DE, and riso, 
are the parameters related to the primordial scalar power spectrum in Eq. ([1]). 

In our calculations, we take the total likelihood to be the products of the separate likelihoods (A) of CMB, LSS 
and SNIa. Defining Xl i = ~21og>Ci, we then have 



2 _ 2 

XL,total — Xl,CMB 



2 I 2 

Xl,LSS I XL,SNIa 



(5) 



If the likelihood function is Gaussian, Xl coincides with the usual definition of up to an additive constant corre- 
sponding to the logarithm of the normalization factor of C. 



^ Available at; |http:/ /cosmologist.info/cosmomc^ 
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The data used for current constraints include the five-year WMAP (WMAP5)[i| as well as some small-scale CMB 
measurements, such as CBI VSAiSj], BOOMERanCfs^ and the newly released ACBAR data^]. For the Large 
Scale Structure information, we use the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample [34|. 
The supernovae data we use are the recently released "Union" compilation of 307 samples 0]. In the calculation 
of the likelihood from SNIa we marginalize over the relevant nuisance parameter [35| . Furthermore, we make use 
of the Hubble Space Telescope (HST) measurement of the Hubble parameter Hq = lOOh km s^^ Mpc""'^ [s^ by 
multiplying the likelihood by a Gaussian likelihood function centered around h = 0.72 and with a standard deviation 
a = 0.08. We also impose a weak Gaussian prior on the baryon density ili,h^ — 0.022 ± 0.002 (1 a) from the Big Bang 
Nucleosynthesis js^l, and a cosmic age tophat prior as 10 Gyr < to < 20 Gyr. 

For the future data, we consider the measurements of LSS from LAM0ST[8j], the CMB from PLANCK 0] and the 
SN la from 5-year SNLS[il. For more information about the mock data, we refer to |25| . 1381 1 . 

For the simulation of LAMOST, we mainly consider the 31? matter power spectrum of galaxies. The simulated 
error of the matter power spectrum including the statistical errors due to sample variance and shot noise are given 

where V is the survey volume and n is the mean galaxy density. In our simulations, we set the redshift of the LAMOST 
main sample to be z ~ 0.2, and the survey area to be 15000 deg^. The total number of galaxies within the survey 
volume is 10'^[|. The maximum k we consider is ~ 0.1 /i Mpc ^. For the simulation with PLANCK, we follow 
the method given in our previous paper We mocked the CMB TT, EE and TE power spectrum by assuming 
the certain fiducial cosmological model. For the detailed techniques, we refer to our previous paper [25|. We have 



also simulated 500 SN la according to the forecast distribution of the SNLS [40|. For the error, we follow the Ref . |41 1 
which takes the magnitude dispersion 0.15 and the systematic error <Tsys — 0.02 x z/1.7. The whole error for each 
data is given by: 



/ 15^ 

0.n,aaa{Zi) = \ <T%siZi) + , (7) 

where rii is the number of supernova of the i'th redshift bin. We take ACDM model as the fiducial model in simulating 
the data. 



III. NUMERICAL RESULTS 

A. Oscillating Primordial Power Spectrum 

The oscillating behavior of the primordial spectral index will imprint on CMB and LSS via the following two 
equations: 

P„(k)=T\k)V^(k), (8) 

where T{k) and Ti{k) are the transfer functions which denote perturbation evolution of matter and photon from 
reheating era up to now [13, H^. Based on the parametrization in Eq. ([T]), in Fig. [Tl we display the matter power 
spectrum Pm{k) and the TT power spectrum C;. The parameters are chosen to be [uso, namp, 7t/re] = [0.96, 0.10, 3.0], 
which can fit the current observational data at 2(t C.L.. The blue solid lines are given by the specific oscillatory power 
index, while the ordinary smooth parametrization are given by the red dashed lines. One can find the modulated 
oscillating structure in the whole linear regime of the matter power spectrum, and this scenario as well as the related 
theoretical model has been considered in Ref.[l3]. Unlike Trans-Plankian mechanism or step- like models which just 
give a local feature in the k space, this kind of signature has effects in the whole linear regime, however it might not 
be easy to distinguish from baryon acoustic oscillation in small scales. From Figll]one can see the effects on CMB 
TT power spectrum with the change of the amplitude of the Doppler peaks. In short, these distinct hints on the 
observations will feed back constraints on the parameters. 

In Table [H we give the main results of our numerical calculation. We present the 1, 2a constraints from the current 
data which include WMAP5, SDSS-LRG matter power spectrum as well as the "union" SN la data. We also provide 
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FIG. 1: Primordial power spectrum (top left), matter power spectrum (top right) and CMB TT power spectrum(bottom) of the 
smooth parametrization (i.e. parameterized by a power law with Us = 0.96 and as — 0) (red) and oscillating parametrization 
(blue). The oscillating parameters are [riso, namp, n/re] = [0.96, 0.10, 3.0] which are within lo. 
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0.974 ± 0.016 ± 0.033 


0.003 


^amp 


0.0 ±0.066 ±0.116 


0.028 
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< 6.256 


< 1.651 


Wo 


n Qco+ti.()!«+ti.i(ii 


0.040 


Wamp 


U.UOU_o. 130-0.276 


0.140 



TABLE I: Constraints on the cosmological parameters from the current and future observations. Here we show the mean values 
and 1,2(7 error bars. For the future measurements, we give the standard deviation of these parameters. For some parameters 
that are only weakly constrained we quote the 95% upper limit. 

the constraints from the simulated PLANCK TT, TE and EE power spectrum, the LAMOST matter power spectrum 
and the 5-year SNLS data. For the future constraints, we mainly give the standard deviation. 

In Fig [21 we give the 1-D probability distribution of n^o, 'T-amp and rifrc. The black solid line is the constraints from 
the current observational data. We get n^o = 0.974 ± 0.016(1ct), and also the constraints on the amplitude of the 
power index |namp| < 0.116 at 95% confidence level. The scale invariant power spectrum is within 2cr C.L., but the 
constraints on the frequency of power index is still very weak, namely, the 2a upper limit is n/re < 6.256. Therefore, 
the current data still allow the oscillating structures on the primordial scalar spectrum. 

Since the present data do not give very stringent constraints on the parameters, it is worthwhile discussing whether 
future data could determine these parameters conclusively. From the Table HI one can see that the constraints get 
tightened significantly. The standard deviation of riso, riamp are 0.003 and 0.028 respectively, while the 2a upper limit 
of Ufre is 1.651. Since mocked data are generated from the fiducial model with the standard power-law primordial 
spectrum, the best fit values of riamp and n/re are expected to be around zero. With this simulation, we find the 
scale invariant power spectrum can be tested at Scr C.L., and the error of the oscillating amplitude can be shrunk by 
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FIG. 2: One-dimensional constraints on the namp(top panel) and n/^e (bottom panel). The black solid lines are given by the 
current data, while the blue dashdot lines are given by the simulated data:SNLS+PLANCK+LAMOST. 



a factor of 2.35. 



B. Oscillating dark energy 

The main contribution of dark energy on CMB is on the geometrical angular diameter distance to the last scattering 
surface. The oscillating EoS can also modulate the TT power spectrum as well as the matter power spectrum, since 
oscillating EOS of DE can leave somewhat similar imprints as oscillating primordial spectrum. For the evolution of 
EoS given in Eq.Q, when wq deviates significantly from —1 and w is matter-like, the contributions to large scale 
CMB are significant due to the Integrated Sachs- Wolfe (ISW) effects. 

Firstly, with current data, we get wq ~ —0.958 ± 0.098, Wamp = O.OSO^Q'iJgQ at 68% confidence level, however, 
AC DM remains a good fit. 

In Fig. [3] we display the two dimensional cross-correlated constraints and one dimensional probability distribution 
of wq and Wamp- The black solid lines are from the current data sets. For the future simulated data, the wq will be 
constrained tightly, however Wamp gets broaden a little bit, which is due to the degeneracy between the wq and the 
Wamp- Quantitatively, the future simulated data shrinks the error bar of wq by a factor of 2.4 as delineated by the blue 
dashdot lines. And interestingly, it also rotates the contour in the WQ-Wamp plane. To understand this phenomenon we 
notice that physics which is interesting to us is the possible deviation of the w from w = —1. In our paramctrization, 
Wq + Wa gives the largest value for this deviation. In Fig. |4l we plot the 1-D probability distribution of wq + Wa- 
From this figure one can see indeed the future data give a tighter constraint. 

Finally, in figure [5] we present the one dimensional constraint on the evolution of w{z) from the current data. One 
can see from this figure AC DM remains a good fit, however the dynamical models with oscillated feature are not 
excluded. 



IV. SUMMARY 

In this paper we have studied the signature of the modulated primordial scalar power spectrum and the oscillating 
EoS of DE in the observations, such as CMB, LSS and SN la. Based on the parameterizations of ns{k) and w{z) 
in Eq.© and Eq. Q, we firstly present the constraints with the current observations. Furthermore, we consider 
the constraints from the future astronomy surveys, for example, the future CMB project PLANCK, the LAMOST 
telescope and 500 SN from the coming 5 year SNLS. 

Our results show that the current data have put constraints on the featured primordial scalar power spectrum 
scenario, however, the limits are rather weak. From global data analysis, we get Ugo = 0.974 ± 0.016 which indicates 
a red power spectrum consistent with our previous analysis. Within 2(7 C.L., with the current observations we have 
I'T'ampI < 0.116 and njre < 6.256, and for the future observations, the error bar can be shrunk significantly 

On the other hand, for the constraints on DE, we find the oscillating EoS of DE can fit the current data well, as 
well as the AC DM model. The amplitude of the oscillation of EoS are limited to be Iwampl < 0.232 at 95% confidence 
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FIG. 3: One-dimensional constraints on the uiq (bottom right panel) and if amp (top panel). The black solid lines are given by 
current data, the blue dashdot lines are given by mocked data:SNLS+PLANCK+LAMOST. 



0.1 - 




-l.a -1.6 -1.4 -1.2 -1 -0.8 -O.B -0.4 



FIG. 4: One-dimensional constraints on the wq -h t/^amp. Different colored line represents different datasets used as before. 
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FIG. 5: Constraints on w{a) = + Wamp sin[^ from the current observations. Median value (central red dotted line), 

68% (area between black solid lines) intervals are shown. The red dotted line is the cosmological constant boundary. 
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level. The future data will give a stronger constraint, especially on the wq parameter. 
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